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Recently, based on the method of collective variables the statistical field theory for multicomponent inhomo- 
geneous systems was formulated [0. Patsahan, I. Mryglod, J.-M. Caillol, Journal of Physical Studies, 2007, 11, 
133]. In this letter we establish a link between this approach and the classical density functional theory for 
inhomogeneous fluids. 
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The powerful tools for the study of equilibrium and non-equilibrium properties of many-particle in- 
teracting systems are those based on the functional methods. In many cases the partition function of such 
systems can be re-expressed as a functional integral after performing the Hubbard-Stratonovich (HS) 
transformation GlU a simple device proposed in the SOies. Nearly at the same time another method, the 
method of collective variables (CVs), that allows one in an explicit way to derive a functional represen- 
tation for many-particle interacting systems, was developed |3, 4]. The method, proposed initially in the 
1950s Q-Hl for the description of the classical charged many particle systems and developed later for the 
needs of the phase transition theory 1 6—10], was in fact one of the first successful attempts to attack the 
problems of statistical physics using the functional integral representation. The CV method is based on: 
(i) the concept of collective coordinates being appropriate for the physics of the system considered (see, 
for instance, |ll|]) and (ii) the functional integral identity 



exp 



[F[p]) = f&p6^[p- p] exp [F[p]) (1) 



valid for classical systems and permitting to derive an exact functional representation for the configura- 
tional Boltzmann factor. Being applied to the fluids, the CV method uses the idea of the reference system 
(RS), one of the basic ideas in the liquid state theory I I 211 . 

Recently, the rigorous scalar field KSSHE (Kac-Siegert-Stratonovich-Hubbard-Edwards) theory (3 
Hill , which uses the HS transformation, was developed to describe the phase equilibria in simple and 
ionic fluids. As was shown 1 15, 16], both theories (KSSHE and CVs) are in close relation. 

Another valuable theoretical approach, which has been extensively employed to study the structural 
and thermodynamic properties of inhomogeneous systems is the classical density-functional theory (DFT) 
(for an overview we refer to (l7|,[l8|]). The central quantity of DFT is the Helmholtz excess free energy ex- 
pressed as a functional of the single-particle density. There are few systems for which this functional is 
known exactly. Thus, successful application of DFT critically depends on judicious choice of an appropri- 
ate Helmholtz energy functional suitable for the system under investigation 1 12] . The main goal of this 
letter is to constitute a link between the CVs based theory and the classical DFT. 

Let us consider the general case of a classical m-component system consisting of N particles among 
which there exist N\ particles of species 1, N2 particles of species 2, ... and N m particles of species m.The 
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potential energy of the system is assumed to be of the form 

= "E Z UafituTj) + E E Vafo), (2) 

a,;8 1V7 a=l i=\ 

where U a p{rj,rj) denotes the interaction potential of two particles and the second term is the potential 
energy due to external forces. 

The pair interaction potential U a ptri,Tj) can be considered as the sum 

U a p(r { ,rj) = vlp(ri,Tj) + w a p(jc u *i), (3) 

where i/°^(r,-,rj) is a potential of a short-range repulsion which in general describes the mutual impen- 
etrability of the particles, while w a p(ji,rj), on the contrary mainly describes the behaviour at moderate 
and large distances. The system with the interaction potential v° a „(r,-,rj) can be regarded as the reference 
system (RS). The fluid of hard spheres is most frequently used as the RS in the liquid state theory since its 
thermodynamic and structural properties are well known. 

Introducing the microscopic density of the ath species in a given configuration 

N a 

Pair) = Y J S{r-r i ) 

i=i 

we can present the grand canonical partition function of the system as follows Hill : 



S[{v a }] = E TH E ^r- E ttS f (dr)ex P 



RS P ~ 

-p y Ni ..N m -^<Pa\w afi \pp) + (V a \pa) 



(4) 



with (df) = n a ^Fn,,, drjv a = dr"dr" . . . dr^ , being the element of the configurational space of N parti- 
cles. In the right hand side of equation (|4) Dirac's brackets notations 

E / drdr' p a (r)w a p(r,r')pp{r') = (p a \ w a p\pp), 
E / dry/ a {r)p a [r) = (y/ a \p a ) 

a J 

are introduced and summation over repeated indices is meant. In l|4), Y^ s N denotes the contribution 
from a m-component RS, 

v a {r) = v a + vl-py/ a {r) (5) 

is the local chemical potential of the ath species, v a = P\i a - 3 In A a , A a is the de Broglie thermal wave- 
length and is the self-energy of the ath species = pw aa (r,r)/2. For a given volume V, E[{v a }] is a 
function of the temperature T and a log-convex functional of the local chemical potentials v a [r). 

Using Q} we can present the Boltzmann factor which does not include the RS interaction in the form 

expl^iPalWaplPf})^ = J@p@w exp|i(p a |w a/5 |p / g) + i(w a |{p a -p a })j. (6) 

Here p a (r) is the collective variable which describes the field of the number particle density of the ath 
species.The functional integrals which enter the above equation can be given a precise meaning in the 
case where the domain of volume V occupied by particles is a cube of side L with periodic boundary 
conditions which will be implicitly assumed henceforth II 5l Il9fl . This means that we restrict ourselves to 
the fields p a {r) and d) a [r) which can be written as Fourier series. 

Inserting equation (6} in the definition 0) of the grand canonical partition function one obtains an 
exact functional representation 

H[{v a }] = J Sip&o exp{-,Jf>[{v a },{p a ,c»a\]) , (V) 
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where the action J€ [{v a },{p a ,w a }] of the CV field theory reads as 

3€\{v a },{p a ,Ua\\ = -(Palw a j3lp0}-i(«Jp a }-lnE R s[{v a -iw a }] . (8) 
In (8) Hrs [{v a - ioi a }] is the grand canonical partition function of a m-component RS defined as follows: 
S RS [{v*}] = £ -L £ -L. £ J-j f (dr)exp(-^« s ^ + ^1^)) (9) 

with 

v«(r) = v B (r)-ua B (r). (10) 

Some comments are in order. It should be emphasized that the description fTl-flOl is based on the two 
sets of variables {p a \ and {w a } and valid for repulsive, attractive as well as arbitrary pair interactions. 
One can distinguish the following alternative approaches to the application of equations dTt— flOt: 

• Let W denote the matrix of elements w a p (r, r') . If W is the positive-definite matrix, the integration 
over {p a } can be easily performed. As a result, one arrives at the same functional representation as 
that obtained by means of the HS transformation. Thus, we stress that the HS transformation has a 
narrower region of applicability compared to the CV method. 

• Let matrix Wbe arbitrary. In this case, one can start with integration over {d) a }. In general, Hrs[{v*}] 
cannot be calculated exactly. In order to develop a perturbation theory we present the logarithm 
of the grand partition function of the RS in the form of a cumulant expansion 

lnHnsKv*}] = £^J- £ [dl...(dnm ai „ an a,...,ri)<o ai a)...a>a n m. (ID 

n»0 n - ai,...,a„J J 

where i = r, and di = dr,-. In fTO the nth cumulant 9Ji ai ...a n (l, . . . , n) is equal to the rc-particle partial 
truncated (connected) correlation function of RS at v* (r) = v a (r). Substituting (ll) in we obtain 



H[{v a }] - SRs[{v a }] J@p®w expj -~{p a \w a p\pp) 

(-i)" f f 

+i(w a |p a )+£ — — £ / dl... / dnW ai ... an (l,...,n)a) ai m...a) an W 

n»l n - ai,...,a n J J 



The calculation of correlation functions of RS, 9Kai...a„(l>"-> n)> is a separate task. If RS is a hard 
sphere mixture, one can use the fundamental-measure theory fioll and the Percus-Yevick or Carna- 
han-Starling approximations in the non-uniform and uniform cases, respectively. 

Equation Jl2) can be evaluated in a systematic way using the Gaussian distribution as a basic one. 
In particular, using the Gaussian averages one can develop a loop expansion of the grand partition 
function as it was done recently for a one-component fluid fljl . 

• Another way of integrating over {a> a } is to use the steepest descent method. If RS is a mixture of 
ideal gases, the integration in (9} can be performed exactly. We consider this special case in detail. 

Let RS be a m-component mixture of point particles that corresponds, in turn, to the condition Y^ s N =0 
in equation (9). In this case the right hand side of (9) can be easily calculated 

1 -I -I m f I N a 

Hid[{v a }] = £ — £ — ... £ TT-f ]1 / <...dr« a exp £ v a (r«) 



- n 



a=\ 

Using f!3l we can rewrite the action l[8) as follows: 



£ Tnff drexp(v*(r)) =exp[£( drexp(v;(r)) 

N a »0 ^a- \J I J I a J 



(13) 



^[{Va},{Pa,ti)a}\ = ~ (pal W a p\Pp) ~ K^alPa) ~ £ / drexp (v a (r) - iw a (r)J . (14) 
* a J 
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In (14) we take into account the equation fLO) . It is worth noting that in a uniform case, the action dl4l 
coincides with the corresponding expression obtained in 1 23] . 

In order to integrate in fTJ over the CV fields w a [t) we use the steepest descent method. In this case 
the equation 



= 



leads to the relation 



Pair)- exp [v a (r) - iw a (r)) . 
Introducing the notation p a = exp(-ln + v s a ) we can rewrite Sl5\ as 

- SB a ir) = In - ftp a + /3y/ a (r). 

Pa 

Substituting (16) in equation (14) we finally get for S [{v a }] 

H[{v a }] = J Sip exp{-^[{v a },{p a }]) , 
where the action J€ [{p a }] has the form: 



(15) 



(16) 



(17) 



m{v a },{pa}]=Y. [ dr PaW 
a J 



, Pair) 
In— 

Pa 



+ ~ {Pa\Waf}\Pp) ~ {Pa\PPa} + (p a \Pfa)- (18) 



It is remarkable that the contribution from quadratic fluctuations of A<x> a = a> a - c5 a vanishes in the ther- 
modynamic limit 1 24]. Therefore, the functional representation of the grand partition functional given 
by fTD— (18) is exact for the model with the interaction potential (2)-l[3) under condition v^ir) = 0. This 
contradicts the assumption made in (H] about a purely phenomenological character of functional Hamil- 
tonians which are employed in the field-theoretical formalism. 

For a one-component case, equations dl7)-(IT8) coincide with the functional integral representation 
obtained in ElEaliTll . 

Based on (18) one can formulate the MF theory. To this end, from the stationary-point condition 



dJe[{v a },{ Pa }] 



dp a 







Pa = pa 



we obtain the following set of equations for the MF density 

lnp a (r) = v a ir) - Pj^ \ dr' w a pir,r')ppir'). 
P J 

Then, from (17)-fT8). the MF grand potential reads 

f P 
In H M f [{v a }] = £ I dr pair) + - (p a \w a p\pp) . 

a J ^ 

The Helmholtz free energy of the system defined as the Legendre transform 

F[{pa}] = (pa\Pa) ~ jS" 1 InH [{v a }] , 

has the following form in the MF approximation 



(19) 



^MF [{Pa}] = P l Y,\ dr l°« 
a J 



(r) 



In 



Pair) 

Ha 



■-(Pa\w a p\Pp) + {PaWa) ■ 



Introducing an intrinsic free energy & flBl by the relation & = F - (p a \tya) one can present the MF 
intrinsic free energy of the system as follows: 



^MF [{Pa}] = j6 -1 E f dr P« 
a J 



(r) 



In 



Pair) 

Ha 



+ -(Pa\Waf>\Pf))- 



(20) 
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Note that &mf [iPa]] is the functional of the local density p a (r) and has no explicit dependence on the 
external potential y/ a (r). Functional (20\ has the form of the approximate free energy functional used in 
the DFT when the excess free energy functional arising from the interactions is treated at the MF level. 
In this case the variational principle 

lr —\&[{pa}]-(pa\{Va-Va))}=0 
8 Pa(T) 1 J 

leads to the equation (19) and, in turn, to a coincidence of the MF density p a (r) with the equilibrium local 
density. 

Summarizing, we have derived the exact field theoretical representation (equations (17)— (18)) for 
the m-component inhomogeneous system that does not include the hard sphere interaction. Within the 
framework of the MF formulation of the theory we have found the functional of the free energy which is 
analogous to that used in the MF DFT. 

It is worth noting that the long-range interactions within the framework of the DFT are usually treated 
within the MF approximation. By contrast, the CVs based theory enables one to develop the perturbation 
scheme in order to take into account the correlation effects. Moreover, the latter theory can be applied to 
the systems with more complicated interactions than those considered in this work. 

We have also demonstrated that the density functional integral formulation derived in Iz3,l26l - l28ll can 
be obtained from the exact CVs functional representation in some special case. 

We note that in general the integration over {co a } can be performed exactly in thermodynamic limit 
if the RS describes noninteracting particles. When short-range interactions are included into the RS the 
method of steepest descent will produce additional contributions describing the correlations between 
particles in the RS. In particular, this is the case of a hard sphere RS that will be considered elsewhere. 
However, even for RSs of noninteracting particles there are several models where the method devel- 
oped above can be very useful. One of such examples is connected with the statistical field theory of 
anisotropic fluids where the RS may be considered as the system of noninteracting rigid rotators (see, for 
instance, El ). 



References 

1. Stratonovich R.L., Sov. Phys. Solid State, 1958, 2, 1824. 

2. Hubbard J., Phys. Rev. Lett., 1959, 2, 77; doi |10.1103 /PhysRevLe tlT77l 

3. Zubarev D.N., Dokl. Acad. Nauk SSSR, 1954, 95, 757 (in Russian). 

4. Yukhnovsky I.R., Zh. Eksp. Ter. Fiz., 1958, 34, 379 (in Russian). 

5. Bohm D., Pines D., Phys. Rev., 1951, 82, 625; doi 10.1103/PhysRev.82.625 

6. Yukhnovskii I.R., Phase Transitions of the Second Order: Collective Variables Method. World Scientific, Singa- 
pore, 1987. 

7. Yukhnovs'kii I.R., Rivista del Nuovo Cimento, 1989, 12, No. 1, 1; doi |10.1007/BF027 40597 

8. Yukhnovskii I.R., Proceedings of the Steklov Institute of Mathematics, 1992, 2, 223. 

9. Patsagan O.V., Yukhnovskii I.R., Teor. Mat. Fiz., 1990, 83, 72 (in Russian). 

10. Yukhnovskii I.R., Patsahan O.V., J. Stat. Phys., 1995, 81, 647; doi |10.1007/BF0217925l| 

11. Yukhnovskii I. R., Holovko M.F., Statistical Theory of Classical Equilibrium Systems. Naukova Dumka, Kiev, 1980 
(in Russian). 

12. Hansen J.P., McDonald I.R., Theory of simple liquids. Academic Press, 1986. 

13. CaillolJ.-M., Mol. Phys., 2003, 101, 1617; doi 10.1080/0026897031000068488 

14. Caillol J.-M., J. Stat. Phys., 2004, 115, 1461; doi 10.1023/B:JOSS.0000028066.25728.cf 

15. Caillol J.-M., Patsahan O., Mryglod I., Physica A, 2006, 368, 326; doi 10.1016/j.physa.2005.11.010 

16. Patsahan O., Mryglod I., Condens. Matter Phys., 2006, 9, 659. 

17. Evans R., Adv. Phys., 1979, 28, No. 2, 143; doi 10.1080/00018737900101365 

18. Singh Y, Phys. Rep., 1991, 207, 351; doi |10.1016/0370- 1573(91)9009T6l 

19. Patsahan O., Mryglod I., Caillol J.-M., Journal of Physical Studies, 2007, 11, No. 2, 133. 

20. Rosenfeld Y, Phys. Rev. Lett., 1989, 63, 980; doi 10.1103/PhysRevLett.63.980 

21. Lebowitz J.L., Phys. Rev., 1964, 133, 895; doi |l0Ti 03/PhysRev.l33.A895 

22. Mansoori G.A., Carnahan N.F., Starling K.E., Leland T.W., J. Chem. Phys., 1971, 54, 1523; doi 10.1063/1.1675048 



24001-5 



0. Patsahan, I. Mryglod 



23. Yukhnovskii I.R., Honopolskii O.L., Preprint of the Institute for Theoretical Physics, ITP-74-93P, Kiev, 1974, (in 
Russian). 

24. Frusawa H., Hayakawa R., Phys. Rev. E, 1999, 60, R5048; doi 10.1103/PhysRevE.60.R5048 

25. Evans R., Mol. Phys., 1981, 42, 1169; doi 10.1080/00268978100100881 

26. Woo H.-J., Song X., J. Chem. Phys., 2001, 114, 5637; doi 10.1063/1.1353553 

27. di Caprio D., Badiali J.P., J. Phys. A: Math. Theor., 2008, 41, 125401; doi 10.1088/17 51-8113/41/1271 25401 

28. di Caprio D., Holovko M.F., Badiali J.P., Condens. Matter Phys., 2003, 6, 693. 

29. Holovko M., di Caprio D., Kravtsiv I., Condens. Matter Phys., 2011, 14, 33605; doi |10.5488/CMP.14 .33605 



MeTOfl KOTieKTi/iBHUx 3MiHHi/ix: 3B'5i30K 3 Teopieto <J>yHKij;ioHa/iy 
rycTMHi/i 

0. riai4araH, I. Mpnnnofl 

IhtCTHTyT c(>i3MKW KOHfleHCOBaHHX CUCTeM HAH YKpaiHM, By/1. CBEHL4iL4bKOrO, 1, 7901 1 M. /lbBiB, YKpaiHa 

HeflaBHO, BHKopncTOByK)4H MeTOfl KO/ieKTHBhtux 3MiHHHx, 6y/io C(()opMy^bOBaHO CTaTWCTHKO-no^bOBy TeopiK) 
Pflft 6araT0K0Mn0HeHTHHx HeoflHopiflHwx cwcreM [0. Patsahan, I. Mryglod, J.-M. Caillol, Journal of Physical Studi- 
es, 2007, 11, 133]. B u,bOMy noBiflOM/iem-ii mm BcraHOB/iioeMO 3b'a30k Mix \iy\w niflxoflOM i K/iacn4H0K) Teopiro 
cfiyHKujOHa/iy rycTWHW fl/ia HeoflHopiflHux n/ini-iiB. 

K/iK>HOBi c/iOBa: <pyHKU,iOHajibHi Meiop,vi CTamcnwHo'f cf>i3M<n, Merofl KO/ieKmBHtix 3mihhhx, leop'tn 
<pyHKU,iOHiny rycTi/iHW, 6araroKOMnoHeHTHa HeoflHopiflHa cmcieMa 
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